The reversibility of the Goos-Hanchen shift near the band-crossing structure of 
one-dimemsional photonic crystals containing left-handed metamaterials 
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We perform a theoretical investigation on the Goos-Hachen (GH) shift in one-dimensional photonic 
crystals (IDPCs) containing left-handed metamaterials (LEIMs). We find an unusal effect of the GH 
shift near the photonic band-crossing structure, which is located at the condition, — fci^'dA = 
= TTiTT (m = 1,2,3- ■■), under the inclined incident angle, here A denotes the LHM layer 
and B denotes the dielectric layer. Above the frequency of the band-crossing point (BCP), the GH 
shift changes from negative to positive as the incident angle increases, while the GH shift changes 
reversely below the BCP frequency. This effect is explained in terms of the phase property of the 
band-crossing structure. 
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It is well known that the Goos-Hachen (GH) shift refers 
to the lateral shift of a well-coUimated light beam to- 
tally reflected from the interface of two different media 
[U, 13] ■ Recently the GH shifts have extensively studied 
in various situations, for example, dielectric slab systems 
d, Hi, left-handed materials (LHMs) Otto config- 

uration 0, one-dimensional photonic crystals (IDPCs) 
metal surfaces IDPCs with a nonlinear meta- 
materials [llj , coherent driving atomic systems , and 
electro-optic crystals [l^ and so on. Meanwhile, Li et al. 
have discovered that multilayered structures containing 
LHMs possess zero averaged refractive index gap (zero- 
n gap) [3], which is quite different from a traditional 
Bragg gap. Shadrivov et al. [l5'| have further found that 
in such structures there exists unusual angular depen- 
dencies of the beam transmission in which two photonic 
bands may touch under a certain condition. In this re- 
port, we consider the GH shifts near the band-crossing 
structure of the IDPCs containing the LHMs, and show 
that the GH shifts have unusual properties: Above the 
frequency of the band-crossing point (BCP), the GH shift 
changes from negative to positive as the incident angle 
increases, however the GH shift has the opposite behavior 
below the BCP frequency. 

For the simplicity, we only consider a TE-polarized 
light beam incident from vacuum into the IDPC at an 
incident angle 0, as shown in Fig. 1. The results are sim- 
ilar for TM-polarized light fields with frequencies near 
the band-crossing structure. The IDPC is consisted of 
alternative layers A and B, with the structure {AB)^ , 
where N is period. cLa and ds are widths of layers A and 
B, respectively. The layer A is a lossless LHM with the 
effective dielectric permittivity and magnetic permeabil- 
ity 
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where ojep/2Tr = 10 GHz and Wmp/STr = 10.5 GHz. The 
layer B, without loss of generality, is vacuum [es = 1 
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FIG. 1: Schematic of a [AB)^ structures. Ar and At are the 
reflected and transmitted GH shift. 



and = 1) in this brief report. In general, the electric 
and magnetic fields at any two positions z and z + lS.z in 
the same layer can be related via a transfer matrix [Tg} 



Afj-(Az,cj, ( 



where kl 



_ ( cos [fc^Az] sin [fc^Az] 

iqj sin [/c^Az] cos [fc^Az] 
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of wave vector kj in the jth layer {j 

and c is the light speed in vacuum. Then 



is the z component 
A,B), Qj = 
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the reflection and transmission coefficients [r{(jj,0) and 
t(Lj,9)] can be readily obtained from the transfer matrix 
method [lal 



t{LU,6) 
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where go = cos 9 for the vacuum of the space z < before 
the incident end and the space z > L after the exit end 
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{L is the total length of the IDPC), and Xij = 1, 2) 
are the matrix elements of Xn{uj,9) — JJ Mj{dj,uj,9) 



which represents the total transfer matrix for the finite 
IDPC. For obtaining the GH shift of the incident beam 
with a sufhciently large beam waist (i. e., the beam with 
a very narrow angular spectrum, Ak << k), the GH 
shifts of both the reflected and transmitted beams can 
be expressed explicitly as [1] 
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where 4>^ ^ are the phases of the reflection and transmis- 
sion coefficients, 6r = r{L0,9) and Qt = t{L0,9), and A 
correspond to the wavelength of the incident beam with 
angular frequency uj. 

In order to know the information of the photonic band 
gap for an infinite periodic structure (N —> oo), according 
to Bloch's theorem, the dispersion at any incident angle 
follows the relation ^] 



cos [P^D] = cos 



ki'^^dA + ki^^B 
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where (3^ is the z component of Bloch wave vector, and 

D — dA + dB- The condition of Eq. © having no real 
solution for jS^ is |cos[/3^D]| > 1, which is well-known as 
the Bragg condition of the photonic band gap. In the 
LHMs (i. e., layers A) ki'^^ is negative, while in layers B 
kz is positive. Therefore Eq. ^ becomes 



cos[/3^Z?] — cos 
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When -ki^^dA = ki^Ue ^ rmr {m = 1,2,3,- ••), Eq. 
([7]) has no real solution. Therefore there exists a dis- 
tinctive band gap, the so-called zero-averaged refractive 
index gap (zero-Ti gap) [14]. However, when 

- ki^^^A = ki^^B = mn (to = 1,2,3, •••), (8) 

the upper and lower photonic bands may touch together. 
This phe nomena has been pointed by Shadrivov et al. in 
Ref. [l5|. Here we would like to emphasize that, in the 
case of normal incidence, Eq. ^ leads to the touch ef- 
fect of both the upper and lower photonic bands, however 
such an effect does not belong to the band-crossing effect 
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FIG. 2: (Color online). Photonic band structures of infinite 
(AB)^ structures: (a) dA = ds = 2.5cm, (b) dA = ds = 
7cm, and (c) dA = ds ~ 1.2cm. Heavy solid areas correspond 
to the pass bands, and white areas are the photonic band 
gaps. 



because the upper and lower bands are still parabolic 
shape. In the case of inclined incidence, Eq. ([5]) leads 
to the band-crossing effect, which refers to that both the 
upper and lower photonic bands touch together and the 
dispersion near the touch point is linear. In the follow- 
ing discussion, we focus ourselves to investigate the GH 
shift near the band-crossing structure when the IDPC's 
structure satisfies the condition of Eq. ([5]) in the case of 
inclined incidence. 

In Fig. 2, we plot the photonic band structure in the 
parameter plane (w, fcy), where the transverse wave num- 
ber ky = ^ sin 9 is relative to incident angle 9 for a fixed 
u!. From Fig. 2(a) and (b), it is seen that the band- 
crossing effect occurs at the angles where Eq. ^ is sat- 
isfied for the inclined incidences. Close to the BCP, the 
dispersion relation between uj and (ky — kym) is linear, 
where kym is the transverse wave number for the mth 
band-crossing structure. We find that such band-crossing 
structures are similar to the Dirac band structure in some 
two-dimensional photonic crystals, which is consisted of 
triangular or honeycomb lattices of the dielectric cylin- 
ders ti4[l9|,|20j,l2lll2^. Recently, we also found that it is 
possible to realize such a band-crossing structure in ho- 
mogenous negative-zero-positive index media [1^. Near 
the band-crossing structure, one has demonstrated some 
unusual properties, such as conical diffraction (19j . static 
diffusive property of the light fields [2l|, [2^ and Zitterbe- 
wegung of optical pulses [2^ [23|. For comparison, when 

—ki'^^dA = ki^^dB < TT, we plot Fig. 2(c) to demon- 
strate a zero-rl gap, which is an omnidirectional gap and 
is almost independence of incident angle {25^ . 

In Fig. 3(a), we plot the typical GH shifts of trans- 
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FIG. 3: (Color online), (a) The typical GH shift A as a 
function of angle 6 near the band-crossing structure with m = 
1. (b) Contour plot of (a) as functions of frequency and angle. 
The GH shifts near the band edges of (c) the lower passing 
band and (d) the upper passing band at both sides of the 
zero-n gap in Fig. 2(c). 



mitted light beams near a band-crossing structure with 
m = 1, as shown in Fig. 2(a). The results for the re- 
flected light beams are similar and are not plotted in all 
following cases. Fig. 3(a) shows that, above the BCP 
frequency {uj > lobcp, in our cases ujbcp ~ 27rx 7.24137 
GHz), there is a negative GH shift as incident angle 9 
transits from the gap region into the passing band region, 
however it becomes a positive GH shift as 9 increases 
from the passing band into the gap region. Therefore the 
GH shift for the light with frequencies above the BCP 
frequency can change from negative to positive with the 
increasing of 9. Conversely, below the BCP frequency 
(w < ijJbcp), the GH shift changes from positive to neg- 
ative with the increasing of 9. Figure 3(b) is a contour 
plot of Fig. 3(a). It is clear that, near the band-crossing 
structure, the GH shift can be changed from positive to 
negative or from negative to positive, depending on light 
frequency located below or above the BCP frequency. For 
comparison, we plot Figs. 3(c) and 3(d) to show the GH 
shifts in the case of Fig. 2(c). It is seen that the GH shifts 
near the band edges of the zero-n gap are always posi- 
tive, which is similar to that in the conventional IDPCs 
but is very different from the present cases. Therefore 
we conclude that the GH shifts inside the band-crossing 
structures are very unusual: it can change from negative 
to positive above the BCP frequency and change con- 
versely below that frequency. We can expect that there 
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FIG. 4: (Color online). The GH shift A near the band- 
crossing structures with (a) m = 3 and (b) m = 2 as in 
Fig. 2(b). 



are similar GH effects near the band-crossing structures 
with m — 3 and m — 2, see Fig. 4(a) and 4(b), which 
correspond to the cases in Fig. 2(b). From Fig. 3(a), 
and Figs. 4(a) and 4(b), it is further found that the GH 
shift changes almost linearly from positive to negative 
(or from negative to positive) inside the passing bands of 
such band-crossing structures. 

To gain a deeper insight into the physical mechanism 
for such reversed GH shifts near the band-crossing struc- 
tures, we plot Fig. 5 to show the relative phase change 
as a function of 9 at different frequencies. From Fig. 5, 
the relative phase change of the transmitted light beam, 
respect to the phase of the incident light beam, varies 
from negative to positive as the light frequency crosses 
over the BCP frequency. Here we have to point out that 
the BCP frequency corresponds to the frequency with a 
zero-averaged refractive index of the total IDPC. The rel- 
ative phase is always positive above the BCP frequency; 
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FIG. 5: (Color online). The relative phase shift as a function 
of angle Q at both sides of the band-crossing structure with 
m = \. On each curve, the range between two solid squared 
points denotes the range of the passing band. 
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meanwhile it varies initially from a small to large positive 
value when 6 moves from the gap region into the passing 
band and then deeroascs from the maximal phase change 
to a small one when moves out of the passing band 
again. Thus the GH shifts for the light beams with fre- 
quencies above the BCP frequency change from negative 
to positive, and the GH shifts can be very large negative 
(or large positive) near the transition regions from the 
gap region (or the passing band) to the passing band (or 
the gap region). From Fig. 5, it is clear that below the 
BCP frequency, the relative phase is completely opposite 
with that above the BCP frequency. Therefore the GH 
shift is reversed near the band-crossing structure when 
the light frequency is located below the BCP frequency. 
Finally, it should be pointed out that the phase change 
with respect to is a quadratic-like curve inside the pass- 
ing band (see the regions between two squared points on 
each curve in Fig. 5), therefore the shifts inside the pass- 
ing bands change linearly from negative to positive above 



the BCP frequency (or from positive to negative below 
the BCP frequency). 

In summary, we have investigated the GH shifts in 
IDPCs containing LHMs. It is found that there is an 
unusual effect for the GH shift near the photonic band- 
crossing structure, which appears at the resonant condi- 
tion —ki'^^d.A = ki^^ds = TTin (m = 1,2,3- ••) at the 
inclined incident angle. Above the BCP frequency, the 
GH shift changes from negative to positive as 9 increases, 
while the GH shift changes reversely below the BCP fre- 
quency. Such unusual GH shifts are helpful for under- 
standing on the special properties of the band-crossing 
structures. Finally we have explained this effect in terms 
of the phase properties of the upper and lower passing 
band of the band-crossing structure. 

Acknowledgments: This work is supported by 
NSFC05-06/01, CUHK 401806 and 2060360, and FRG 
of HKBU, and NSFC (10604047). 



[1] F. Goos and H. Hachen, Ann. Phys. 1, 333 (1947); F. 

Goos and H. Hanchcn, Ann. Phys. 5, 251 (1949). 
[2] K. Artmann, Ann. Phys. 2, 87 (1948). 
[3] C. F. Li, Phys. Rev. Lett. 91, 133903 (2003). 
[4] L. G. Wang, H. Chen, and S. Y. Zhu, Opt. Lett. 30, 2936 

(2005). 

[5] P. R. Barman, Phys. Rev. E 66, 067603 (2002). 

[6] A. Lakhtakia, Electromagnetics 23, 71 (2003). 

[7] I. V. Shadrivov, A. A. Zharov, and Y. S. Kivshar, AppL 

Phys. Lett. 83, 2713 (2003). 
[8] D. Felbacq, A. Moreau and R. Smaai, Opt. Lett. 28, 1633 

(2003). 

[9] L. G. Wang and S. Y. Zhu, Opt. Lett. 31, 101 (2006). 
[10] P. T. Leung, C. W. Chen, and H. P. Chiang, Opt. Com- 

mun. 276, 206 (2007); 
[11] R. R. Wei, X. Chen, J. W. Tao, C. F. Li, Phys. Lett. A 

372, 6797 (2008). 
[12] L. G. Wang, M. Ikram, and M. S. Zubairy, Phys. Rev. A 

77, 023811 (2008). 
[13] X. Chen, M. Sheng, Z. F. Zhang, and C. F. Li, J. Appl. 

Phys. 104, 123101 (2008). 
[14] J. Li, L. Zhou, C. T. Chan, and P. Sheng, Phys. Rev. 

Lett. 90, 083901 (2003). 
[15] I. V. Shadrivov, A. A. Sukhorukov, and Y. S. Kivshar, 



Appl. Phys. Lett. 82, 3820 (2003). 
[16] L. G. Wang, N. H. Liu, Q. Lin and S. Y. Zhu, Phys. Rev. 

E,70 016601 (2004). 
[17] M. Contini, C. Sibilia, M. Scalora, G. D'Aguanno, M. 

Bertolotti, M. J. Bloemer, C. M. Bowden, and I. Nefedov, 

Phys. Rev. E 60, 4891 (1999). 
[18] M. Plihal and A. A. Maradudin, Phys. Rev. B 44, 8565 

(1991). 

[19] O. Peleg, G. Bartal, B. Froodman, O. Manela, M. Scgcv, 
and D. N. Christodouhdos, Phys. Rev. Lett. 98, 103901 
(2007). 

[20] F. D. M. Haldane and S. Raghu, Phys. Rev. Lett. 100, 
013904 (2008). 

[21] R. A. Scpkhanov, Y. B. Bazaliy, and C. W. J. Beenakker, 

Phys. Rev. A 75, 063813 (2007). 
[22] X. Zhang, Phys. Rev. Lett. 100, 113903 (2008). 
[23] L. G. Wang, Z. G. Wang, J. X. Zhang, and S. Y. Zhu, 

Opt. Lett. 34, 1510 (2009). 
[24] L. G. Wang, Z. G. Wang, and S. Y. Zhu, Europhysics 

Lett, (accepted) (2009). 
[25] H. T. Jiang, H. Chen, H. Q. Li, Y. W. Zhang and S. Y. 

Zhu, Appl. Phys. Lett. 83 5386 (2003). 



